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NUMERICAL  SIMULATION  OF  A  LOW  DENSITY 
PLASMA  EROSION  OPENING  SWITCH 

I,  Introduction 

The  presence  of  plasma  in  the  vacuum  output  transmission  line  of  a  pulsed 
power  generator  can  enhance  or  degrade  the  power  pulse  delivered  to  the  vacuum 
load.  The  uncontrolled  generation  of  high  density  plasma  in  the  vacuum  line 
generally  degrades  the  coupling  of  the  generator  to  the  load  by  providing  an 
alternative  path  for  the  generator  current.  In  the  worst  case,  the  entire 
generator  current  can  be  shunted  away  from  the  1  oad  by  the  pi  asma.  Sources  of 
uncontrolled  plasma  include  intense  UV  irradiation  of  surfaces  and  plasmas 
from  explosive  field  emission.  However,  by  injection  of  low  density  plasmas 
into  the  vacuum  line,  the  alternative  current  path  to  ground  through  the 
plasma  can  be  controlled  as  a  function  of  time,  and  the  power  pulse  delivered 
to  the  load  intentionally  modified. 

In  particular,  the  geometry  of  the  vacuum  transmission  line  and  the 
properties  of  the  Injected  plasma  can  be  chosen  so  that  a  plasma  erosion 
opening  switch  (PEGS)  configuration^  is  obtained,  as  illustrated  in  Fig.  1. 
Originally  this  configuration  was  successfully  used  as  a  plasma  erosion  switch 
(PES)  for  prepulse  suppression, ^  which  suggested  its  potential  as  an  opening 

-I 

switch. As  an  opening  switch  the  PEOS  can  be  used  for  inductive  storage 
applications.  The  switch  plasma  conducts  the  generator  current  for  some 
period  of  time,  while  the  energy  of  the  generator  is  stored  inductively  in  the 
vacuum  transmission  line  between  the  generator  and  the  plasma.  At  a 
predetermined  time,  the  plasma  will  cease  to  act  as  a  1  ow  impedance  path  to 

ground,  and  the  current  will  be  swiftly  diverted  to  the  load.  This  inductive 
Manuacript  approvad  Novambar  4, 1986. 


storage/PEOS  technique  modifies  the  output  of  the  generator  to  deliver  a 
similar  current,  higher  voltage,  shorter  duration  pulse  to  the  load.  The 
generator  can  charge  the  vacuum  inductor  at  a  voltage  less  than  the  load 

voltage,  reducing  voltage  stress  on  the  generator  vacuum  insulator  and 
limiting  the  uncontrolled  production  of  plasma  in  the  vacuum  transmission 

line.  High  voltage  at  the  load  exists  only  after  the  PEGS  opens  at  a  time 
when  most  of  the  generator  energy  is  already  stored  in  the  vacuum  section. 
Significant  power  and  voltage  multiplication  (compared  with  the  case  of  no 
PEGS)  is  obtained  only  if  the  opening  time  of  the  PEGS  is  much  shorter  than 
the  conduction  phase  of  the  PEGS.  The  proper  plasma  conditions^  are  crucial 
for  achieving  the  required  fast  opening  time. 

The  work  presented  in  this  paper  treats  the  conduction  of  current  through 
12  -3 

a  low  density  (  IG  cm  )  collisionless  plasma  injected  between  two  coaxial 
conducting  cylinders  as  illustrated  in  Fig.  1.  The  plasma  is  injected  through 
the  anode  towards  the  cathode  with  a  flow  velocity,  Vp,  and  is  assumed 
azimuthally  symmetric.  The  opening  process,  when  current  is  diverted  to  a 
load,  is  also  treated.  The  2  1/2  D  electromagnetic  particle-in-cell  code, 
MASK,^  is  used  to  simulate  the  problem.  Current  is  driven  through  the  plasma 
so  that  the  IGG  kA  level  is  reached  in  5  ns.  Space  charge  limited  electron 
emission  is  allowed  on  the  conducting  surfaces  and  particles  encountering  the 
surfaces  are  absorbed,  A  partially  transmitting  boundary  condition  at  the 
load  end  is  used  to  simulate  a  terminating  load  impedance. 

Results  show  that  initially  a  narrow  radial  current  channel  forms  at  the 
generator  side  of  the  plasma.  The  current  carrying  electrons  include  both 
plasma  electrons  and  electrons  emitted  from  the  cathode.  The  electrons  are 
bent  slightly  toward  the  load  by  the  ^  force  as  they  travel  toward  the 
anode  where  they  are  absorbed.  Generally  their  early  trajectories  are 
predominately  radial,  but  the  slight  bending  of  their  orbits  produces  an  axial 


displacement  of  electrons  (and  space  charge  buildup)  toward  the  load.  This 

results  in  a  layer  of  negative  space  charge  on  the  load  side  of  the  current 

channel  and  a  layer  of  positive  space  charge  on  the  generator  side  of  the 
current  channel  and  thus  a  strong  axial  electric  field  in  the  current 

channel.  As  the  current  increases,  electrons  can  then  _Ex^  drift  radially 

across  the  magnetic  field  lines.  The  axial  electric  field  is  shorted  out  at 

the  anode  and  cathode. 

At  the  cathode  the  induced  radial  electric  field  is  enhanced  by  the 
electric  field  associated  with  the  space  charge  region  in  the  current 

channel.  This  large  radial  field  produces  the  necessary  electron  emission 
from  the  cathode  and  radially  accelerates  these  electrons  into  the  current 
channel.  At  the  anode  the  induced  radial  electric  field  is  reduced  by  the 
electric  field  associated  with  the  space  charge  region  in  the  current 

channel .  The  anode  end  of  the  current  channel  propagates  downstream  toward 
the  1 oad,  whil  e  the  cathode  end  of  the  channel  where  el  ectrons  enter  the 
channel  (emission  site),  moves  only  slightly  (depending  only  on  the  £  x  ^ 
force).  Thus  the  channel  orientation  begins  by  being  radial,  but  becomes  more 
and  more  diagonal  as  the  current  increases. 

The  migration  of  the  anode  end  of  the  current  channel  is  caused  by 
electrons  x  drifting  axially  toward  the  load  and  neutralizing  the  ion 
space  charge.  These  electrons  enter  the  system  from  the  anode  and  fill  in 
behind  the  current  channel  on  the  generator  side  as  it  moves  toward  the 
load.  When  the  current  channel  reaches  the  load  end  of  the  plasma,  the  axial 
space  charge  electric  field  can  no  longer  be  maintained.  At  this  time  the 
electron  flow  is  quickly  magnetically  insulated  and  the  current  is  diverted  to 
the  load.  While  the  plasma  conducts  current,  the  current  channel  remains 
narrow  (i.e.,  much  less  than  the  plasma  width)  and  the  current  density 


increases.  The  plasma  flow  velocity  has  little  effect  on  the  time  of  opening 
or  current  density.  The  width  of  the  channel  is  determined  by  the  amount  of 
space  charge  buildup  required  to  produce  the  proper  electric  field  and  the 
associated  £xB_ drift  velocity  for  current  conduction.  In  the  case  of  low  mass 
density  (e.g.,  when  the  ion  species  is  protons),  significant  ^x8_  axi al 
displacement  of  the  emission  site  and  current  channel  also  occurs. 

As  set  up,  this  problem  is  designed  to  be  the  first  step  in  simulating  a 
low  density  lO^^cm"^),  high  flow  velocity  (-^  lO^cm/sec)  PEGS.  Scaling 
from  present  experimental  results®  to  high  voltage  applications^  of  the  PEGS, 
such  as  for  PBFA  II,  leads  to  consideration  of  PEGS  operation  in  this  low 
density,  high  flow  velocity  regime.  However,  because  the  simulations  treat 
this  new  regime,  the  physical  processes  controlling  the  PEGS  operation  may  not 
be  the  same  as  observed  in  the  higher  density  experiments.  In  fact,  the 
results  presented  here  and  in  previously  reported  work,®  are  not  consistent 
with  either  the  experimental  results  at  higher  density  h  10  cm  )  or  the 
theoretical  model ^  which  has  been  used  to  explain  these  higher  density  PEGS 
experiments.  For  example,  experimental  results^*®  show  that  the  current 
channel  broadens  as  the  current  increases  while  the  current  density  remains 
approximately  constant.  The  switch  opens  when  the  current  channel  width 
reaches  the  full  width  of  the  plasma.  The  current  density  in  the  plasma  is 
proportional  to  the  ion  flux  density,  n^Vp,  at  the  cathode.  These  results  are 
not  reproduced  in  the  simulations.  There  are,  however,  fundamental 
differences  in  the  two  problems.  In  the  simulations:  a)  the  density  is  an 
order  of  magnitude  lower;  b)  the  dl/dt  is  a  factor  of  five  larger  (which  is 
more  similar  to  the  current  risetime  of  PBFA  II);  c)  the  timescale  is  an  order 
of  magnitude  shorter;  d)  the  plasma  electron  temperature  is  more  than  two 
orders  of  magnitude  higher  and  e)  azimuthal  symmetry  is  imposed.  The  presence 


of  anomalous  collisions  due  to  instabilities  could  reconcile  the  experimental 
and  simulation  results  by  allowing  electrons  to  cross  magnetic  field  lines  and 
carry  current  radially  across  the  switch.  This  would  eliminate  the  need  for 
the  large  space  charge  separation  and  electric  fields  found  to  be  required  in 
these  simulations  for  current  flow.  In  addition  the  emission  would 

spread  further  across  the  cathode  than  in  the  present  simulations  because  the 
local  enhancement  of  emission  due  to  the  large  space  charge  separation  would 
not  be  present.  This  collisional  current  carrying  mechanism  could  give 
results  more  similar  to  those  observed  experimentally;  however,  it  is 
precluded  from  appearing  in  the  present  simulations  because  of  the  differences 
a)  through  e)  discussed  above.  This  mechanism  will  be  studied  in  future  work. 

In  Sec.  II  the  numerical  simulations  are  discussed  in  detail,  including 
all  assumptions  and  boundary  conditions.  In  Sec.  Ill  the  results  of  the 
simulations  are  presented.  Results  are  shown  for  different  ion  species, 
different  densities,  different  plasma  injection  velocities  and  different 
current  risetimes,  A  detailed  discussion  of  the  results  is  given  in  Sec.  IV 
and  Sec.  V  includes  a  discussion  of  the  conclusions  that  can  be  drawn  from 
this  work.  Plans  for  extending  this  work  to  include  the  effects  of  anomalous 
collisions  are  also  discussed  in  the  last  section. 

II.  Numerical  Simulations 

The  problem  to  be  simulated  deals  with  the  conduction  of  current  through 
a  low  density  (  ^  10  cm  )  collisionless  plasma  injected  between  two  coaxial 
conducting  cylinders  as  shown  in  Fig.  1,  The  opening  process,  during  which 
current  is  diverted  to  the  load,  is  also  dealt  with.  This  is  the  first  step 
in  simulating  the  low  density  operation  of  a  PEOS.  The  2  1/2-0 

electromagnetic  particle-in-cell  code  MASK  is  used  to  properly  treat  both  the 


switch  plasma  and  the  boundary  emission  of  electrons  into  the  plasma.  The 
problem  is  done  in  r-z  geometry  and  azimuthal  symmetry  is  imposed. 
Simulations  were  also  run  in  planar  geometry,  ^side  from  geometfry  effects, 
the  current  conduction  and  opening  process  were  similar;  therefore  only  the 
cylindrical  geometry  cases  will  be  discussed  here. 

The  coaxial  electrodes  and  the  switch  plasma  form  part  of  a  circuit  with 
a  generator  at  the  left  and  a  load  in  parallel  with  the  switch  at  the  right. 
The  generator  in  the  simulations  is  modeled  as  a  current  source  with  current 
rising  linearly  from  0  to  100  kA  in  about  5  ns.  The  anode  and  cathode  are 
perfect  conductors,  while  the  boundary  at  the  load  is  given  an  "impedance"  by 
specifying  a  transmission  coefficient  for  electromagnetic  waves  incident  on 
the  radial  surface.  The  cathode  radius  is  2.5  cm,  and  the  anode  radius  is  5 
cm.  The  simulation  region  is  12.8  cm  long  axially,  with  the  plasma  occupying 
the  volume  between  z  =  2.5  cm  and  z  =  7.5  cm.  The  regions  on  both  sides  of 
the  plasma  between  the  plasma  and  the  axial  boundaries  are  vacuum  regions. 
The  cathode  and  anode  boundaries  are  also  the  sources  of  ions  and  electrons. 
Plasma  electrons  and  ions  are  injected  down  from  the  anode  between  z  =  2.5  cm 

Q 

and  7.5  cm  with  flow  velocity  Vp  =  1.8x10°  cm/s.  Electrons  are  also  allowed 
to  be  emitted  from  both  the  cathode  and  anode  as  required  by  space  charge 
limited  emission.  Particles  striking  any  of  the  boundaries  are  absorbed.  The 
initial  plasma  fill  itself  is  also  given  a  fl  ow  velocity  toward  the  cathode  of 
Vp  =  1.8x10®  cm/s  and  has  a  density  of  0^  =  0^  ~  ''e  ~  cm"®.  A  cold  start 

is  used  whereby  cold  electrons  and  ions  are  initially  placed  directly  on  top 
of  each  other  and  given  a  radially  inward  flow  velocity. 

The  simulation  region  is  divided  into  a  128  (in  z)  by  50  (in  r)  mesh, 
which  is  uniform  in  each  direction.  The  mesh  sizes  are  then  Ar  =  0.05  cm  and 
Az  =  0.1  cm.  Four  macroparticles  of  each  specie  are  used  initially  in  each 
cell  in  the  plasma  region.  Simulations  run  with  different  numbers  of 


macroparticles  per  cell  showed  that  this  was  about  the  minimum  number  needed 
to  maintain  good  particle  statistics.  Larger  numbers  of  macroparticles  per 
cell  produced  similar  results.  The  skin  depth  for  the  plasma  with  n^  * 
lO^^cm"^  is  about  c/Upg  =  0.5  cm,  and  is  easily  resolved  in  all  the 
simulations  to  ensure  accurate  attenuation  of  the  fields  by  the  plasma.  In 
addition,  the  simulation  region  is  broad  enough  to  prevent  the  fields  from 
leaking  through  to  the  load  before  opening. 

Some  numerical  questions  of  concern  include  the  resolution  of  the  Debye 

length  and  the  correct  treatment  of  the  electron  emission  region.  The 

experimentally  measured  value  for  the  initial  plasma  temperature  in  a  PEOS 

configuration  is  about  Tp  =  5  eV.®  This  would  require  a  prohibitively  small 

mesh  size  to  resolve  the  Debye  length.  Since  it  is  important  to  simulate  the 

full  switch  region  in  order  to  investigate  both  current  conduction  and  switch 

opening,  an  artificially  high  plasma  temperature  is  used.  This  temperature  is 

specified  by  equating  the  mesh  size  with  the  Debye  length  for  a  given  plasma 

density.  Using  a  cold  start  this  temperature  is  quickly  reached  (<  1  ns)  by 

numerical  heating,  which  ends  after  the  Debye  length  increases  to  the  mesh 

size.  Simulations  performed  with  an  initially  warm  plasma,  Tg 

2  2 

=  Aitn^e  (Ar)  ,  show  identical  mechanisms  of  current  conduction  and  switch 
opening.  This  test  verifies  that  the  cold  start  can  be  used  since  the  early 
numerical  heating  does  not  influence  the  general  switch  behavior. 

The  artificially  high  plasma  temperature  and  corresponding  large  mesh 
size,  however,  can  have  an  important  effect  on  the  simulation  results.  The 
large  mesh  size  eliminates  smaller  wavelength  instabilities  and  waves,  and  the 
high  temperature  can  reduce  the  growth  rate  and  even  stabilize  those  streaming 
instabilities  with  wavelengths  that  can  be  resolved.  The  imposition  of 

azimuthal  symmetry  could  also  suppress  certain  instabilities  because  only 
m  =  0  modes  are  allowed.  These  initial  simulations  are  not  intended  to 
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address  the  effect  of  plasma  instabilities  on  switch  behavior;  however,  a 
comparison  of  the  simulation  results  with  experimental  data  demonstrates  the 
importance  of  further  investigations  of  the  role  of  plasma  instabilities  and 
anomolous  resistivity.  This  will  be  left  for  future  work. 

The  large  mesh  size  used  in  the  simulations  also  sets  a  limit  on  the  size 

of  the  space  charge  limited  emission  gap,  d,  that  can  be  resolved  at  the 

conducting  boundaries.  For  a  given  gap  voltage,  this  in  turn  limits  the 

3/2  2 

available  emitted  current  density  to  J  =  a  4*^  /d  .  The  mesh  size  in  the  gap 
normal  direction  is  chosen  in  the  simulations  to  resolve  the  Debye  length 
associated  with  the  energy  of  emitted  electrons.  Thus  the  emission  of 
electrons  is  not  affected  by  numerical  effects.  A  simulation  with  a  smaller 
radial  mesh  size  was  run  to  test  the  numerical  treatment  of  the  emission 
sites.  Results  verify  that  the  emission  processes  are  being  treated  accurately. 

In  the  next  section,  the  results  of  the  simulations  will  be  discussed. 
Simulations  have  been  run  at  densities  of  Oq  =  SxlO^^cm"^  and  at  IxlO^^cm"^. 
The  current  risetime,  tp,  has  also  been  varied  with  current  reaching  100  kA  in 
5  ns  and  10  ns.  Results  using  both  H'*'  and  C'*’  ions  are  discussed  as  well  as 
results  using  different  plasma  flow  velocities.  Parameters  for  the  runs  that 
will  be  discussed  are  given  in  Table  I,  An  analysis  of  switch  behavior  for 
each  run  will  be  given  and  compared  with  the  other  runs. 


Table  I:  Parameters  for  each  of  five  simulation  runs. 


Run  # 
1 
2 

3 

4 

5 


^5 

n  (cm"  ) 

1x10^2 

0.5x10^2 

1x10^2 

1x10^2 

1x10^2 


tp  (ns) 
5 
5 

10 

5 

5 


Speci es 


Vp  (cm/us) 
180 
180 
180 
180 
0 


Ill,  Results 

This  section  will  be  devoted  to  the  presentation  of  results  obtained  from 
the  computer  simulations.  Simulations  of  the  five  different  cases  shown  in 
Table  I  are  discussed  here,  but  the  physical  interpretation  of  the  results  is 
left  for  the  next  section.  The  different  simulations  have  many  aspects  in 
common.  Switch  operation  is  qualitatively  the  same  for  all  five  cases,  with 
differences  appearing  in  such  things  as  the  timing  of  switch  opening  and  the 
speed  of  current  channel  motion  downstream.  Four  identifiable  aspects  of 
switch  operation  are:  1)  early  current  channel  formation,  2)  movement  of  the 
anode  end  of  the  channel  toward  the  load,  3)  movement  of  the  cathode  end  of 
the  channel  toward  the  load,  and  4)  switch  opening  by  magnetic  insulation  of 
the  electrons.  These  four  processes  can  be  seen  in  the  series  of  Figs.  2-4. 
These  plots  are  snapshots  taken  of  various  computational  diagnostics  at  times 
t  *  0.4,  3.0,  4.0,  5.0  and  6.0  ns  from  Case  1  in  Table  I,  which  will  be  used 
as  the  model  problem  here.  Plots  between  0.4  and  3.0  ns  have  been  omitted 
because  very  little  activity  occurs  there  and  because  the  numerical  heating 
mentioned  in  Section  II  obscures  the  results  between  0.6  and  1.2  ns. 

In  Fig.  2  current  density  vectors  are  plotted  in  the  computational 
region.  At  t  =  0.4  ns  (Fig.  2a)  current  flows  between  the  anode  and  cathode 
at  the  generator  side  of  the  plasma.  The  current  remains  confined  in  a 
channel  that  never  broadens  more  than  about  1  cm  as  it  penetrates  through  the 
plasma.  The  channel  evolves  in  time  as  seen  in  Figs.  2b-2e.  The  anode  end  of 
the  current  channel  propagates  towards  the  load  through  the  plasma  region 
faster  than  the  cathode  end.  Eventually,  as  shown  in  Fig,  2e  current  no 
longer  crosses  between  anode  and  cathode  in  the  region  of  the  plasma  but  has 
been  diverted  to  the  load.  Switch  opening  is  defined  as  the  point  when  90%  or 
more  of  the  generator  current  has  been  switched  to  the  load.  For  Case  1  in 
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Table  I,  the  switch  conducts  current  for  about  5  ns  before  reaching  this 
opening  point. 

The  evolution  of  the  electrons  and  ions  in  the  switch  is  shown  in  Figs.  3 
and  4.  The  electrons  in  Fig.  3  include  the  original  plasma  electrons,  the 
emitted  electrons  and  the  electrons  injected  from  the  anode  at  the  plasma  flow 
velocity.  Early  in  time  (t  <  0.4  ns),  electrons  at  the  generator  side  of  the 
plasma  are  accelerated  toward  the  anode,  leaving  a  gap  in  the  lower  left 
corner.  This  gap  remains  occupied  by  the  plasma  ions,  whose  space  charge 
elicits  electron  emission  into  the  gap.  The  electrons  at  the  generator  side 
of  the  plasma  at  t  =  0.4  ns  have  been  displaced  slightly  toward  the  load, 
thereby  exposing  a  layer  of  ions  along  that  edge.  At  all  times  the  current 
across  the  switch  is  carried  predominantly  by  the  electrons,  so  that  the 
electron  particle  plots  can  be  seen  to  follow  the  current  channel  plots 
described  earlier.  Current  carrying  electrons  are  emitted  from  the  cathode 
and  either  enter  the  current  channel  directly  (because  they  have  been  emitted 
near  the  cathode  end  of  the  channel )  or  enter  the  current  channel  after 
flowing  downstream  along  the  cathode  in  a  magnetically  insulated  manner  until 
they  reach  the  current  channel.  The  orbits  of  these  electrons  will  be  shown 
later  in  Fig.  5a.  Two  special  features  of  Fig.  3  should  be  noted.  As  the 
anode  end  of  the  channel  moves  downstream,  a  diagonal  depletion  layer  is  seen 
in  the  electron  distribution.  This  depletion  layer  evolves  from  the  layer  of 
exposed  ions  seen  along  edge  of  the  plasma  on  the  generator  side  at  t  =  0.4 
ns.  The  layer  develops  as  the  anode  end  of  the  current  channel  moves  toward 
the  load.  By  t  «  3  ns,  the  layer  extends  diagonally  from  the  original  gap 
region  at  the  cathode  on  the  generator  side  all  the  way  to  the  anode  nearly  at 
the  load  end.  This  depletion  layer  grows  in  time,  with  the  current  channel 
forming  its  boundary  on  the  load  side.  The  ion  space  charge  in  the  depletion 


layer  contributes  to  the  electrostatic  electric  field  necessary  for 
the^xB  drift  of  the  electrons  across  the  switch.  Late  in  time  (t  >  5  ns) 
the  electron  motion  is  almost  completely  magnetically  insulated,  with  more 
than  90%  of  the  current  flowing  along  the  cathode  to  the  load. 

The  other  notable  feature  of  the  electron  particle  plots  is  the  block  of 
electrons  in  the  upper  left  hand  corner  that  defines  the  boundary  on  the 
generator  side  of  the  depletion  layer.  These  electrons  consist  largely  of 
electrons  emitted  and  injected  from  the  anode  surface  which  attempt  to 
neutralize  the  ion  space  charge  in  the  depletion  layer.  The  orbits  of  these 
electrons  will  be  shown  later  in  Fig.  5b.  Their  role  in  controlling  the 
motion  of  the  anode  end  of  the  current  channel  will  be  discussed  in  Section 
IV. 

The  electron  behavior  should  be  compared  with  the  ion  particle  plots  seen 
in  Fig.  4.  At  t  *  0.4  ns,  ions  have  not  had  enough  time  to  react  and  are  seen 
in  their  original  configuration.  At  later  times,  a  comparison  of  the  electron 
and  ion  plots  shows  that  although  ions  occupying  the  electron  depletion  layer 
have  a  local  density  much  larger  than  the  local  electron  density,  they  have  a 
reduced  density  compared  with  the  original  plasma  fill.  Also  seen  in  the  ion 
plots  is  the  formation  and  motion  of  a  high  density  layer  of  ions.  This 
structure  is  formed  when  ions  are  accelerated  toward  the  load  by  the  electric 
field  in  the  current  channel  that  produces  the  £xB^  drift  of  the  electrons. 
This  is  a  manifestiation  of  ^x^  accel eration.  This  layer  of  ions  acts  like  a 
leaky  piston,  picking  up  a  certain  fraction  of  ions  as  it  moves  downstream. 
The  ions  left  behind  by  the  piston  occupy  the  electron  depletion  layer 
mentioned  earlier.  The  position  of  the  piston  marks  the  cathode  end  of  the 
current  channel  because  the  ion  space  charge  in  the  piston  stimulates  the 
emission  of  current  carrying  electrons  directly  into  the  channel.  Other 
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electrons  generally  enter  the  channel  after  being  emitted  on  the  generator 
side  of  the  current  channel  and  flowing  along  the  cathode  until  reaching  the 
radial  current  channel.  At  t  =  5  ns,  most  of  the  ions  in  the  piston  have  been 
absorbed  at  the  cathode,  leaving  only  a  small  clump  of  ions  at  the  lower  right 
corner  of  the  plasma.  These  ions  are  responsible  for  the  spray  of  electrons 
in  the  corresponding  electron  plot.  One  final  feature  to  note  is  that  ions 
which  are  continually  injected  at  the  anode  fill  in  behind  the  depletion  layer 
on  the  generator  side  at  the  original  plasma  density. 

Typical  electron  trajectories  in  the  current  channel  are  shown  in  Fig.  5 
at  t  4  ns.  Note  that  the  cathode  end  of  the  current  channel  is  at  4.5  cm 
while  the  anode  end  is  at  7.25  cm.  Figure  5a  illustrates  that  the 
trajectories  of  many  of  the  current  carrying  electrons  are  straight  and  that 
even  those  electrons  which  do  execute  perpendicular  gyromotion  can  cross  the 
A-K  gap  in  very  few  gyroperiods.  The  importance  of  this  observation  is  made 
clear  in  Section  IV.  The  electrons  in  Fig.  5b  are  executing  orbits  in  the 
region  of  plasma  on  the  generator  side  of  the  current  channel.  These 
electrons  attempt  to  neutralize  the  space  charge  in  the  electron  depletion 
layer.  They  E_  x  drift  axially  into  the  depletion  layer,  thereby  reducing 
the  local  ion  space  charge.  Consequently,  the  electron  current  channel  moves 
toward  the  load  in  order  to  maintain  the  space  charge  separation  required  for 
the  E,  X  ^  drift  current  conduction. 

The  differences  in  the  behavior  of  the  switch  for  the  various  cases 
listed  in  Table  I  are  predominantly  ones  of  timing.  All  four  aspects  of 
switch  behavior  discussed  earlier  are  present  in  all  five  cases.  In  Figs.  6 
and  7,  the  relative  motion  of  the  anode  end  of  the  current  channel  and  the 
piston  motion  along  the  cathode  (or  the  motion  of  the  cathode  end  of  the 
current  channel)  are  shown  for  all  five  cases.  Figure  6  demonstrates  that  the 


motion  of  the  anode  end  of  the  channel  Is  approximately  logarithmic  in  time. 
Figure  7  demonstrates  a  nearly  parabolic  motion  for  the  ion  piston.  Analytic 
interpretations  and  explanations  of  these  characteristics  will  be  given  in  the 
next  section. 

Finally,  a  comparison  is  made  in  Table  II  of  the  approximate  times  of 
opening  for  the  different  cases  in  Table  I.  The  time  of  opening  is  defined  as 
the  point  when  90%  or  more  current  is  diverted  to  the  load.  This  can  be 
compared  with  the  time  when  only  about  10%  of  the  current  is  diverted  to  the 
load,  also  shown  in  Table  II.  In  all  respects,  other  than  timings,  the 
current  carried  by  the  switch  plasmas  for  the  various  cases  evolves  very  much 
like  Case  1.  In  Fig.  8,  magnetic  "probe"  data  for  Case  1  is  exhibited  for  two 
probes  at  r  »  4.6  cm.  One  probe  (z  *  2.8  cm)  is  just  inside  the  upstream 
boundary  of  the  original  plasma  configuration,  the  other  lies  close  to  the 
load  (z  *  12  cm).  The  probe  near  the  load  shows  that  the  switch  plasma 
conducts  more  than  90%  of  the  current  presented  to  it  until  about  t  «  2  ns,  at 
which  point  the  switch  can  be  said  to  begin  to  open  and  allow  more  and  more 
current  to  flow  to  the  load.  At  t  *  5  ns  90%  of  the  current  is  flowing 
through  the  load. 


Table  II:  Time  at  which  10%  and  90%  of  the  total  current  begins  to  flow  in 
the  load  for  each  of  the  five  runs  listed  in  Table  I. 


Case 

1 

2 

3 

4 

5 


10%  point  (ns) 
2.0 
1.6 
3.6 
2.8 
2.2 
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90%  point  (ns) 
5.0 

3.5 

6.5 
5.9 
5.1 


In  the  next  section,  the  data  presented  here  is  interpreted  and  analyzed 
to  understand  the  underlying  physics  of  the  switch  behavior  discussed  in  this 
section. 


IV.  Discussion 

In  Sec.  Ill  the  results  of  several  runs  demonstrated  current  conduction 
and  opening  under  different  plasma  conditions.  The  features  that  these 
simulations  display  are  now  discussed  and  compared.  In  Sec.  Ill,  Case  1  in 
Table  I  was  used  as  a  baseline  for  comparison  because  it  best  illustrates  all 
the  features  of  the  switch  behavior.  Although  Case  1  is  still  used  for 
comparisons,  in  this  section  Case  4  in  Table  I  is  used  as  a  baseline  for  much 
of  the  analysis.  Because  of  the  higher  mass  C'*’  ions  in  Case  4,  the  analysis 
of  the  other  physical  effects  can  be  carried  out  without  being  obscured 

by  £  X  ion  motion.  Early  time  current  conduction  and  current  channel 
formation  are  discussed  first.  The  mechanisms  involved  in  current  conduction 
in  this  channel  are  discussed  next,  together  with  an  examination  of  the  role 
that  electron  orbits  play  in  determining  the  width  of  the  channel.  Following 
this  the  evolution  of  the  channel  in  time  is  described 

including  £  x  B_  effects.  This  section  concludes  with  a  comparison  of  the  five 
runs  within  the  context  of  this  discussion. 

Very  early  in  time,  current  is  carried  at  the  generator  side  of  the 
plasma  by  plasma  electrons  accelerated  radially  towards  the  anode  by  the 
induced  radial  electric  field.  The  plasma  prevents  the  fields  from 
penetrating  more  than  about  a  skin  depth  (c/“pg  0,5  cm)  into  the  switch 
region.  The  current  carrying  electrons  are  also  displaced  slightly  toward  the 
load  by  the  ^x^  force,  resulting  in  a  thin  layer  (  <  0.5  cm)  of  space  charge 
and  a  self-consistent  axial  electric  field  as  illustrated  in  Fig.  9a.  The 


massive  ions  are  unable  to  respond  on  a  fast  timescale;  however  on  a  longer 
timescale,  ion  acceleration  in  the  electric  field  is  a  manifestation  of  the 
single  fluid  ^  x  force.  On  the  shorter  timescale  the  axial  electric  field 
provides  an  £  x  8_  drift  channel  on  the  load  side  of  this  space  charge  layer 
for  radial  electron  current  conduction.  The  ion  gyroperiod  is  long  compared 
with  the  time  scale  of  interest  so  that  ion  £  x  ^  drift  motion  is  not 
important.  Electron  trajectories  are  basically  straight,  carrying  the 
electrons  radially  from  the  cathode  toward  the  anode  in  this  channel. 

Because  the  axial  electric  field  must  vanish  at  the  conducting  cathode, 
another  mechanism  must  be  responsible  for  transporting  the  electrons  from  the 
cathode  surface  to  this  drift  channel.  The  initial  radial  flow  of  plasma 
electrons  leave  behind  a  region  of  exposed  ions  or  gap  between  the  plasma  and 
the  cathode  at  the  generator  side  of  the  plasma  as  was  seen  in  Fig.  3a.  The 
space  charge  associated  with  these  ions  locally  enhances  the  Induced  radial 
electric  field  and  hence  the  emission  of  electrons  from  the  cathode.  These 
electrons  are  accelerated  across  this  gap  and  enter  the  £  x  £  drift  channel  to 
be  carried  across  to  the  anode.  In  order  to  cross  the  gap  the  effective 
electron  gyroradius  must  be  large  compared  with  the  gap  size. 

At  early  times  when  electrons  enter  the  £  x  £  drift  channel,  they  flow  in 
a  nearly  force-free  fashion.  That  is  to  say  the  force  from  the  electric  field 
is  nearly  balanced  by  the  ^  x  £  force  so  that  the  £  x  £  drift  velocity  of  the 
electrons  is  much  larger  than  the  velocity  associated  with  oscillations  about 
the  force-free  trajectory  (later  in  time  larger  oscillations  are  observed). 
This  force-free  flow  is  stable  in  a  fluid  sense,  in  that  perturbations  in  the 
current  channel  position  produce  forces  which  restore  the  channel  to  its 
original  position.  If  the  electron  current  channel  were  displaced  toward  the 
load  more  space  charge  would  be  produced  and  the  electric  field  would  be 


increased.  This  unbalanced  electric  field  would  pull  the  channel  back  toward 
the  generator  and  its  original  position.  Similarly  if  the  electron  current 
channel  were  displaced  toward  the  generator  the  space  charge  would  be  reduced 
as  well  as  the  electric  field.  The  unbalanced  _v  x  ^  force  would  now  draw  the 
current  channel  back  toward  the  load  and  its  original  position.  Thus  current 
conduction  provided  by  the  E.  x  8^  radial  drift  of  electron  across  the  plasma  is 
stabl  e. 

Consider  now  the  space  charge  layer  and  x  dri ft  channel  in  more 
detail.  These  structures  are  illustrated  schematically  in  Fig.  9a.  The  width 
of  the  space  charge  layer  is  controlled  by  the  electric  field  required  to 
carry  the  instantaneous  current,  I(t),  The  width  can  be  estimated  by  assuming 
that  an  annulus  of  electrons  of  density  an  is  displaced  toward  the  load  a 
distance  equal  to  its  width,  w.  Here  a(t)  <  1,  The  total  width  of  the 
current  channel  and  the  ion  layer  is  2w,  The  resulting  electric  field  then 
provides  the  electron  drift,  Vq  Ec/8,  for  current  conduction  with 

I  »  2iirwJ  »  2irrwn(a+l)eEc/8.  (1) 

Poisson's  Equation  can  be  used  to  estimate  E  yielding 

E  “  4irnaez  =  E^z/w,  0  <  z  <  w.  (2) 

Here  z  »  0  at  the  load  end  of  the  current  channel  and  increases  moving  toward 
the  generator.  Similarly  Ampere's  Law  can  be  used  to  calculate  8  as 

8  "  Zlzl^rc  =  8^(r)z/w,  0  <  z  <  w.  (3) 
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The  width  of  the  current  channel  is  then 

w(t)  »  I(t)/(2ffrnec[a(o+l)]^^^).  (4) 

The  drift  velocity  can  also  be  expressed  as  =  c(a/l+a)^^^. 

If  a  approaches  unity,  then  Vq  approaches  c//?.  From  the  simulation  with  C"'’ 
(Case  4  in  Table  I)  w  quickly  grows  to  <  0.5  cm  in  the  first  nanosecond  and 

then  grows  slowly  while  o(t)  starts  out  small  and  approaches  unity  as  I(t) 

increases.  At  4  ns,  I  *  80  kA,  and  a  0.35,  so  that  from  Eq.  (2)  w  1.0  cm 

for  n  ®  10^^  cm"^  and  r  =  3.5  cm.  Thus  at  late  times  the  electrons  drift  at 
nearly  c//?.  This  is  in  reasonable  agreement  with  the  simulation  results. 
When  X  ^  acceleration  causes  significant  ion  motion,  n  actually  increases  in 
the  current  channel.  Increased  density  in  turn  reduces  the  required  w 

and/or  a.  This  is  evident  in  the  simulations  at  late  times  (t  -  5  ns)  for 

carbon  ions  in  Case  4,  and  plays  an  important  part  even  at  early  times  for  the 
proton  cases. 

Electron  motion  in  the  E^  x  8^  dri ft  channel  can  be  analyzed 
nonrel  ativistically  by  assuming  static  fields  and  using  a  slab  geometry.  The 
assumption  that  the  fields  do  not  vary  in  time  is  reasonable  because  the 
electron  gyroperiod,  is  ^  0.1  ns  and  because  the  electrons  in  the 

current  channel  cross  from  the  cathode  to  the  anode  in  <  1  ns,  both  of  which 
are  short  compared  with  the  5-10  ns  current  risetimes  in  the  simulations.  The 
nonrel ativi Stic  treatment  and  the  use  of  slab  geometry  are  for  simplicity. 
The  equations  of  motion  for  the  electrons  are  then 


dv,  eE  z  ev  B  z 

—  »  -2-  -  X  °  (fi\ 

dt  m  w  m  cw  *  '  ' 

e  e 

where  Eq  and  are  defined  as  in  Eqs.  (2)  and  (3)  and  the  same  coordinate 
system  is  used  as  in  Eqs,  (2)  and  (3)  with  x  measuring  in  the  direction  from 
the  cathode  to  the  anode.  By  transforming  to  a  frame  drifting  at 
the  ^  X  ^  drift  velocity,  Vg  =  E^c/Bq,  the  electron  motion  is  found  to  involve 
a  constant  drift  velocity  and  a  gyromotion,  v  =  ^  -  Vge^ ,  described  by 


u  V  z 
CO  gz 
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(7) 


d  V 
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where  =  eB^/m^c,  Multiplying  Eq.  (7)  by  Vgj^  and  Eq.  (8)  by  Vg^,  adding 
these  two  equations  and  integrating  yields 
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(9) 


where  Vg^^g  and  Vg^g  are  the  initial  electron  velocity  components.  Integrating 
Eq.  (7)  also  yields 


(d 


gx  gxo 


+  f 7^  7 

"^0 


(10) 


where  Zg  is  the  initial  axial  position  of  the  electron  when  it  enters 
the  E^  X  ^  current  channel.  Combining  Eqs.  (9)  and  (10)  results  in 


gz  dt 
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This  equation  can  be  solved  in  terms  of  elliptic  integrals  and  inverted  to 
yield  z(t).  Here,  however,  the  intent  is  only  to  show  that  self-consistent 
orbits  exist  in  the  current  channel  and  to  compare  the  characteristic  of  these 
orbits  with  those  seen  in  the  simulation  (see  Fig,  5a). 

Electrons  enter  the  current  channel  at  Zq  with  energy 

^  eE^z  eE 

(12) 


,  .  “  eE  z  eE^z: 

Im  v^  >  /  r— ^]dz  =  °° 

Z  eo'^''o'-  w-* 


2w 


where  the  inequality  is  used  because  energy  acquired  from  the  induced  electric 
field  has  been  neglected.  In  general,  energy  acquired  from  the  induced  fields 
is  small  compared  with  the  energy  associated  with  the  electrostatic  potential 
well  of  the  current  channel  in  the  plasma.  The  integral  can  be  evaluated  from 
0  to  Zq  regardless  of  the  actual  path  the  electron  takes  because  the 
electrostatic  force  in  conservative  and  the  cathode  is  at  zero  potential. 

If  an  electron  enters  the  current  channel  with  Vq»Vj^q«Vq  (i.e., 
Vgjjo®Vg2o’01  at  Zg^ZQ,  then  it  will  drift  straight  from  the  cathode  end  to  the 
anode  end  of  the  current  channel.  From  Eq.  (12)  one  find  that 

Zq/w  <  (vq/wu^^ )^2  ^  0.43  for  Vq  ~  l,5xl0^®cm/s,  w  1  cm  and  SxlO^^s"^ 

at  t  »  4  ns  and  r  =  3,5  cm.  Thus  electrons  injected  near  the  center 
of  £  X  drift  current  channel  and  parallel  to  it  will  have  straight  line 
orbits  as  seen  in  Fig.  5a.  Electrons  injected  at  an  angle  or  at  some  other  Zq 
will  gyrate  about  Zq  as  it  travels  from  the  cathode  to  the  anode, 

A 

Electrons  injected  into  the  channel  at  Zg  with  ''q  -^o 

confined  in  the  channel  if  z(t)  is  bounded  by  w  >z(t)>0.  The  turning  points 

or  bounds  of  the  orbit  can  be  found  from  Eq,  (11)  by  setting  v-_  *  0,  For 

2 

example,  if  Vgj^g  ■  0,  the  turning  points  are  at  z  ^  *  (z^  ±  ^wv^^^/w^^)  . 

2  ^2 

Similarly,  if  Vg^g  *  0.  z+  »  Zg  and  z.  »  (z^  -  ^''g^g^^gc^  *  electron 


entering  the  channel  at  Zg  =  w  with  Vg^g  =  0  and  Vg^^g  *  Vg  -  Vq  >  0,  then  has 
turning  points  at  z+  =  w  and  z_  =  [w^  -  4w(Vg-VQ)/u»ggi^2  .  Such  an  electron  is 

o  o  o 

confined  to  the  channel  if  z  >  0  or  v  /vf  <  +  wu  ^/4vp,)  .  Using  the  same 

“  OD  COU"^ 

numbers  as  before  for  Vq,  w  and  this  shows  that  electrons  with  energies 

2 

as  much  as  five  times  that  associated  with  the  drift  motion  80  keV) 

can  be  confined  in  the  channel.  Current  carrying  electron  kinetic  energies  as 
high  as  400  keV  are  observed  in  the  simulations  at  t  =  4  ns. 

The  average  gyroperiod  for  an  electron  in  the  current  channel  is 
T  =  Zirw/u  <z>  ~  4it/u»  ,  where  <z>  is  the  time  averaged  axial  position. 

CO  CO 

Comparing  the  distance  an  electron  drifts,  VfjT,  in  one  gyroperiod  with  the 
axial  excursion,  Az  =  |z^-  z_|,  during  a  gyration,  yields  the  ratio 
VnT/Az  >  4irvr>/w  w  ~  2.3.  Again  the  same  numbers  as  those  used  following  Eq. 
(12)  were  used.  This  result  implies  that  the  orbits  will  be  elongated  in  the 
X  direction  as  observed  in  Fig.  5a  from  the  simulations.  Thus  the  orbits 
found  from  Eqs.  (5)  and  (6)  quite  adequately  describe  the  electron  motion  in 
the  current  channel  observed  in  the  simulations. 

The  electron  orbits  shown  in  Fig.  5b  can  be  used  to  explain  the  motion  of 
the  anode  end  of  the  current  channel.  On  the  generator  side  of  the  switch 
region  electrons  can  enter  the  plasma  from  the  anode  behind  the  ion  layer  and 
the  current  channel.  These  electrons  are  initially  bent  away  from  the  ion 
layer  and  current  channel  toward  the  generator  by  the  magnetic  field. 
However,  once  in  the  plasma,  they  can  £  x  £  drift  axially  toward  the  load  and 
thus  into  the  ion  layer.  Near  the  anode  the£x£drift  is  predominantly 
axial  because  the  E^  field  must  vanish  at  the  conducting  anode.  The  induced 
Ep  then  provides  the  axial  drift.  With  electrons  drifting  into  and 
neutralizing  the  ion  space  charge  layer  from  the  generator  side,  the 
electron  current  channel  must  move  toward  the  load  in  order  to  maintain  the 
space  separation  required  for  the  proper  E^. 
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On  the  generator  side  of  the  current  channel  B  ~  B^t/x  The  induced 
can  be  estimated  using  Faraday's  law  as  B^t/cx,  where  3/3z.  The 
axial  drift  velocity  then  scales  as  E^c/B  ~  Z/t.  Thus  the  position  of  the 
anode  end  of  the  current  channel , 

Z-(t)  ^  /z/t  dt  CL  ln(t),  (13) 

d 

should  have  a  logrithmic  dependence  on  time.  Figure  6  shows  a  plot  of  the 
position  of  the  anode  end  of  the  current  channel  as  a  function  of  time.  The 
behavior  in  time  exhibited  in  the  simulations  is  in  very  good  agreement  with 
this  model.  This  indicates  that  the  anode  end  of  the  current  channel  is 
controlled  by  electrons  drifting  into  the  ion  layer  from  the  generator  side  of 
the  switch  region. 

The  position  of  the  cathode  end  of  the  current  channel  is  controlled  by 
another  mechanism.  As  seen  in  Fig.  4,  a  dense  layer  of  protons  moves  toward 
the  load  and  controls  the  main  emission  site  of  the  electrons  and  hence  the 
current  channel  source.  The  motion  of  this  channel  source  can  be  modeled 
using  a  piston  or  snowplow  approach.  This  structure  acts  like  a  leaky  piston 
because  it  appears  to  sweep  up  most  but  not  all  of  the  ions  in  its  path  and 
carry  them  along  as  it  travels  toward  the  load.  A  simple  analysis  of  this 
structure  in  cylindrical  geometry  can  be  made.  Assume  that  the  total  current 
in  the  channel  of  width  w  is  supplied  by  electrons  of  uniform  density  n(t) 
traveling  with  the  E_  x  B_  drift  provided  by  space  charge  separation.  Then  I(t) 
is  given  by  Eq.  (1)  with  n(a  +  1)  set  equal  to  n(t)  and  r  evaluated  at  the 
cathode.  Ions  are  collected  in  the  piston  and  are  accelerated  by  the  electric 
field  in  the  channel  such  that 
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(14) 


where  M  is  the  ion  mass  and  Xp  the  distance  travelled  by  the  piston.  This  is 
equivalent  to  a  simpl  e  ^  x  ^  accel  eration  of  the  plasma  as  demonstrated  by 
noting  that  enE  =  envpB/c  =  JB/c  where  Vq  is  just  the  el  ectron  E_  x  ^  dri ft 
velocity.  Although  there  is  charge  separation  in  the  piston  (see  Fig.  9b),  it 
is  assumed  that  the  ion  and  electron  densities  averaged  over  the  piston  width 
are  roughly  equal.  Thus,  the  density  of  electrons  and  ions  in  the  channel  is 
only  a  function  of  the  amount  of  ions  swept  up.  It  is  also  assumed  here  that 
only  a  fraction,  B,  of  the  ions  crossed  by  the  piston  are  swept  up.  This  is 
the  leaky  piston  assumption.  Then  the  density  of  current  carrying  electrons 
can  be  written  as 


Bx  (t) 

n(t)  ^  nq(l  -g--),  (15) 

where  n^  is  the  original  plasma  fill  density. 

Once  the  piston  has  begun  to  move,  the  electron  current  channel  and  the 
associated  structures  can  be  pictured  as  in  Fig.  9b.  Plasma  with  the  original 
fill  density,  n^,  exists  on  both  the  load  and  generator  sides  of  these 
structures.  The  current  flows  diagonally  across  the  plasma  in  the  higher 
density  or  piston  region  where  there  exists  negative  space  charge.  On  the 
generator  side  of  this  current  channel  an  equivalent  positive  space  charge  is 
spread  over  the  ion  layer  at  high  density  (n.  >  n^)  and  the  depletion  layer 
(no  >  n^  >>  Og).  The  total  width  of  the  structure  including  the  current 
channel,  the  ion  layer,  and  the  depletion  layer  is  still  ^  2w  as  defined  in 
Eq.  (4). 

Equations  (1),  (14)  and  (15)  can  now  be  combined  to  obtain  an  equation 
describing  the  motion  of  the  piston 
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(16) 


d^x. 


Zif  •  (23F?;‘  6n„M(dXp/dt)  )/*1wn^(l+6Xp/w). 


dt‘ 


Equation  (16)  can  be  solved  given  the  initial  conditions  Xp(0)  and  Xp(0),the 
current,  I(t),  the  magnetic  field  in  the  channel,  B(t),  and  estimates  of  w 
and  6  from  the  simulation.  Ion  acceleration  and  collection  in  the  piston  is 
hampered  if  the  electron  current  channel  sweeps  through  the  plasma  so  fast 
that  ions  cannot  be  picked  up  in  the  "potential  well”  of  the  ^  x  ^  current 
channel.  For  very  fast  sweep  speeds  the  piston  is  very  leaky  (i.e. ,  B  0), 
while  at  slower  sweep  speeds  more  ions  can  be  accelerated.  This  leakage 
occurs  in  all  the  cases  that  were  simulated  because  the  anode  end  of  the 
current  channel  moves  rapidly  through  the  plasma.  The  piston  is  most  leaky  at 
early  times  when  the  anode  end  of  the  current  channel  is  moving  fastest.  The 
axial  acceleration  of  ions  does  not  occur  close  to  the  cathode  because  the 
axial  electric  field  vanishes  there.  Motion  of  the  piston  along  the  cathode 
occurs  because  the  current  channel  extends  diagonally  across  the  plasma  and 
ions  are  accelerated  downward  toward  the  cathode  as  well  as  axially  toward  the 
load.  The  position  of  the  piston  at  the  cathode,  z^(t),  is  then  proportional 
to  its  displacement  Xp(t)  through 


Zc(t)  =  Xp(t)/sin0, 


(17) 


where  9  is  the  angle  made  by  the  piston  with  the  cathode.  Results  of  the 
piston  model  will  be  shown  for  Z(.(t). 

In  order  to  test  the  validity  of  the  piston  model  a  problem  is  chosen 
where  the  motion  of  the  anode  and  cathode  ends  of  the  channel  are  relatively 
decoupled.  The  best  example  is  the  C**  problem  after  the  anode  end  of  the 
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current  channel  has  nearly  reached  the  load  edge  of  the  original  plasma. 
Before  this  point  is  reached,  the  anode  end  of  the  current  channel  moves 
through  the  plasma  producing  little  ion  acceleration  because  of  the  large  ion 
inertia  (see  Figs.  6  and  7).  After  the  anode  end  slows  down  (  4  ns),  ion 

acceleration  can  occur  relatively  unaffected  by  the  motion  of  the  anode  end. 
Equations  (16)  and  (17)  will  be  used  to  model  this  situation  with  6  assumed 
constant.  The  initial  conditions  and  z'^(tQ)  at  t^  =  4  ns  are  taken 

from  Fig.  7.  The  current,  I(t),  magnetic  field  in  the  channel,  B(t),  the 
channel  width,  w,  and  the  "leakage"  factor,  B,  are  all  taken  from  the 
simulation  (Case  4,  Table  I). 

Results  are  shown  in  Fig.  10  for  Zj.(4  ns)  =  3.35  cm,  2'(.(4  ns)  « 
5.5xl0®cm/s,  rQ=2.5  cm,  w=l  cm,  0=0.3,  nQ=10^^cm"^,  and  9=30®.  Figure  10a 
depicts  the  motion  of  the  piston  along  the  cathode,  Z(.(t),  and  shows  good 
agreement  with  the  actual  piston  position  taken  from  the  simulation.  Figures 
10b  and  10c  plot  the  calculated  density  of  the  current  carrying  electrons 
(equal  to  the  density  of  collected  ions)  and  the  transverse  electric  field  in 
the^x^drift  channel.  These  two  quantities  are  more  difficult  to  compare 
with  the  simulation  results  but  there  appears  to  be  good  agreement  as  well. 

The  last  physical  process  to  be  described  involves  the  current  transfer 
to  the  load.  Most  current  carrying  electrons  cross  to  the  anode  in  the 
plasma,  however,  some  electrons  can  leak  out  into  the  vacuum  region. 
Electrons  which  exit  the  current  channel  at  such  an  angle  as  to  carry  them 
into  this  field  free  region  between  the  plasma  and  the  load  will  follow 
ballistic  orbits  until  striking  a  surface.  Ions  are  not  present  in  this 
region  to  provide  the  E^  field  required  for  radial  drifting.  Some  of  these 
electrons  reach  the  load  and  provide  an  initial  load  current.  This  load 
current  then  produces  a  magnetic  field  in  the  region  on  the  load  side  of  the 
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plasma  which  begins  to  magnetically  insulate  more  of  the  electrons  from 
reaching  the  anode  as  they  exit  the  current  channel.  Once  the  anode  end  of 
the  current  channel  reaches  the  load  side  of  the  original  plasma  region,  full 
switch  opening  begins.  Electrons  on  the  generator  side  of  the  current  channel 
continue  to  drift  axially  toward  the  load  and  neutralize  the  ion  space  charge 
in  the  electron  depletion  layer.  Since  there  are  no  further  ions  to  uncover 
on  the  load  side,  the  balance  between  the  £  and  ^x  Enforces  on  the  electron 
current  channel  is  upset  and  the  ^  x  ^  force  dominates.  The  current  channel 
bends  toward  the  load,  and  the  anode  end  of  current  channel  detaches  from  the 
anode.  Electrons  exiting  the  current  channel  are  now  sprayed  into  the  vacuum 
region.  As  the  load  current  increases,  the  electron  current  channel  is  pulled 
toward  the  cathode,  and  ions  in  the  remaining  ion  layer  and  current  channel 
are  also  swept  down  to  the  cathode.  When  the  switch  is  completely  open,  there 
remains:  a  plasma  at  the  anode,  a  gap  with  radial  ion  space  charge  limited 
current  flow,  and  magnetically  insulated  axial  electron  current  flow  along  the 
cathode  surface.  This  final  configuration  is  illustrated  in  Figs.  2e,  3e  and 
4e. 

This  section  closes  with  a  comparison  of  the  results  of  the  five  cases 
from  Table  I.  The  logarithmic  time  dependence  of  the  position  of  the  anode 
end  of  the  current  channel  as  seen  in  Fig.  6  is  expected  from  Eq.  (13). 
Except  for  Case  2,  all  cases  show  roughly  the  same  motion.  Because  of  the 
factor  of  two  lower  density  in  Case  2,  one  would  expect  the  depletion  layer, 
ion  layer  and  current  channel  structure  to  be  broader.  From  Eq.  (4)  this 
width  for  Case  2  should  be  about  twice  as  large  as  that  for  Case  1,  as  is 
confirmed  in  the  simulations.  If  l  in  Eq.  (13)  is  roughly  approximated  by 
this  width,  then  the  velocity  of  the  anode  end  of  the  current  channel  in  Case 
2  should  be  about  twice  as  fast  as  the  other  cases,  in  agreement  with  what  is 
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observed.  The  slightly  faster  notion  of  Cases  1  and  5  over  Cases  3  and  4  is 
most  likely  a  result  of  1  arger  ^  x  ^  piston  motion  in  the  body  of  the 
plasma.  Piston  motion  is  less  in  Case  3  because  of  the  lower  current  and  in 
Case  4  because  of  the  higher  mass. 

The  cathode  end  of  the  current  channel  behaves  as  would  be  expected  from 
the  piston  model  described  earlier.  That  is,  from  Eq.  (16),  the  initial 
acceleration  of  the  piston  is  roughly  proportional  to  the  square  of  the 
current  and  inversely  proportional  to  the  ion  mass  and  density.  It  is 
expected  that  Case  2,  at  half  the  density  of  Case  1,  would  have  about  twice 
the  acceleration.  Similarly,  Case  3  should  experience  about  1/4  the 
acceleration  because  the  risetime  of  the  current  has  been  reduced  by  1/2. 
Case  4  should  show  about  1/12  the  acceleration  because  of  the  higher  mass  of 
the  carbon  ions.  The  differences  between  Case  1  and  Case  5  should  be 
practically  negligible  in  this  context.  All  these  predictions  are  qualitative 
confirmed  in  the  simulations  as  seen  in  Fig.  7.  It  would  be  unrealistic  to 
expect  exact  agreement  between  the  heuristic  arguments  above  and  the  actual 
motion  of  the  piston's  intersection  with  the  cathode  because  of  the  strong 
coupling  of  the  two  ends  of  the  current  channel.  However,  the  agreement  is 
remarkably  good  as  evidenced  by  Fig,  10. 

The  final  comparison  is  the  time  of  opening  for  this  switch.  Table  II 
shows  the  time  at  which  90%  of  the  total  current  has  been  diverted  to  the 
load.  A  comparison  of  these  results  with  the  results  displayed  in  Fig.  6 
shows  a  strong  correspondence  between  the  time  of  opening  and  the  motion  of 
the  anode  end  of  the  current  channel.  Case  2  opens  the  earliest  at  3,5  ns  and 
is  also  the  case  with  the  fastest  motion  at  the  anode.  Similarly  Case  3  opens 
last  and  moves  slowest  at  the  anode.  Thus  it  can  be  concluded  that  the  motion 
of  the  anode  end  of  the  current  channel  determines  the  time  of  opening.  If  it 
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is  assumed  that  t  in  Eq.  (13)  is  given  roughly  hy  w  in  Eq.  (4),  then 
qualitatively  the  switch  will  open  sooner  for  lower  n.  One  also  finds  for  the 
cases  simulated  that  the  time  of  opening  is  basically  independent  of  the 
plasma  flow  velocity,  Vp.  Evidence  for  this  is  that  Cases  1  and  5  are  nearly 
identical.  The  length  of  the  switch,  however,  should  have  a  strong  effect  on 
the  time  of  opening.  The  longer  the  switch,  the  longer  it  will  take  for  the 
anode  end  of  the  current  channel  to  reach  the  load  end  of  the  plasma.  This 
cannot  be  carried  too  far  since  the  motion  of  the  piston  will  eventually 
overtake  the  motion  of  the  anode  end  of  the  current  channel. 


V,  Conclusions 

In  this  paper  we  have  demonstrated  some  of  the  salient  features  of  a 
simulation  of  a  plasma  switch  that  operates  at  low  density  (  -  10  cm  )  and 
in  the  absence  of  any  plasma  instabilities.  The  simulation  is  done  in 
cylindrical  geometry  with  azimuthal  symmetry  and  the  switch  operates  at  a 
plasma  temperature,  which  for  numerical  reasons,  is  near  5  keV,  Current  is 
driven  through  the  plasma  filling  the  region  between  the  anode  and  cathode 
with  electron  current  provided  by  space  charge  limited  emission  of  electrons 
at  the  cathode.  The  initial  plasma  has  a  radial  flow  velocity,  Vp,  toward  the 
cathode  and  plasma  is  continuously  injected  from  the  anode  at  that  flow 
velocity.  There  is  no  evidence  of  plasma  instabilities  or  anomalous 
resistivity  and  the  switch  operates  largely  through  the  mechanisms  of 
electron  £  x  drift  and  space  charge  limited  emission  of  electrons. 

Four  aspects  of  switch  behavior  have  been  identified:  1)  early  current 
conduction  and  current  channel  formation,  2)  motion  of  the  anode  end  of  the 
current  channel,  3)  motion  of  the  cathode  end  of  the  current  channel,  and  4) 
switch  opening.  Analysis  of  the  anode  motion,  cathode  motion  and  electron 


orbits  in  the  channel  has  been  performed,  showing  good  agreement  with 
simulation  results. 

The  current  channel  is  about  1  cm  wide  throughout  the  process.  The 
emission  site  of  electrons  at  the  cathode  varies  in  position  and  width  during 
the  operation  of  the  switch.  Before  the  development  and  motion  of  an  ion 
"piston"  through  the  plasma,  emission  is  concentrated  at  the  generator  side  of 
the  cathode,  with  emission  "enhanced"  by  the  ion  space  charge  exposed  at  the 
generator  side  of  the  plasma.  The  emission  region  widens  slightly  as  the 
channel  bends  diagonally,  but  never  becomes  as  broad  as  the  original  plasma 
region  before  the  piston  has  developed.  This  diagonal  bending  occurs  because 
the  anode  end  of  the  current  channel  moves  toward  the  load  as  electrons  £  x 
drift  axially  into  the  ion  layer  on  the  generator  side  of  the  current  channel. 

When  ions  begin  to  be  accelerated  in  the  channel  and  form  a  piston,  the 
piston  controls  the  load-side  end  of  the  emission  region.  Emission  is  still 
locally  enhanced  by  the  ion  space  charge  at  the  cathode  end  of  the  ion 
layer.  Emission  continues  on  the  generator  side  of  the  piston,  although 
electrons  emitted  there  flow  axially  along  the  cathode  until  they  reach  the 
current  channel , 

Eventually,  the  electrons  "disconnect"  from  the  ions,  i.e,,  the 
trajectories  of  the  electrons  in  the  current  channel  carry  them  into  a  region 
on  the  load  side  of  the  original  plasma  position  where  ions  cannot  provide  the 
axial  field  required  for  radial  £  x^  drifts.  In  this  case  the  electrons  can 
no  longer  cross  the  switch.  At  this  point  current  is  completely  diverted  to 
the  load  and  any  emitted  electrons  become  magnetically  insulated.  The  motion 
of  the  anode  end  of  the  current  channel  determines  the  time  of  opening  which 
is  found  in  the  simulation  to  be  independent  of  the  plasma  flow  velocity.  The 
motion  of  the  anode  end  of  the  current  channel,  and  thus  time  of  opening,  does 


depend  on  density,  moving  faster  for  lower  density.  The  length  of  the  switch, 
however,  also  has  a  strong  effect  on  the  time  of  opening.  The  longer  the 
switch  the  longer  it  will  take  for  the  anode  end  of  the  current  channel  to 


reach  the  load  end  of  the  plasma.  This  scaling  will  be  limited 
by  ^  X  ^  effects  since  the  piston  will  eventually  overtake  the  motion  of  the 
anode  end  of  the  current  channel  In  a  very  long  switch. 

In  future  work,  an  investigation  will  be  made  of  the  presence  of  plasma 
instabilities  and  anomalous  resistivity  In  the  switch  plasma.  Many 
instabilities  cannot  be  simulated  using  the  present  code  because  of  the 
assumption  of  azimuthal  symmetry  and  the  use  of  large  mesh  size  and  high 
plasma  temperature.  However,  it  may  be  possible  to  model  these  instabilities 
by  introducing  appropriate  transport  terms  in  the  particle  equations  of 
motion.  The  effect  of  a  collision  term  has  been  investigated  in  a  steady 
state  fluid  treatment  of  the  problem.  These  preliminary  results  show  that 
switch  behavior  could  be  considerably  altered  under  these  conditions,  giving 
wider  current  channels  and  more  predominantly  radial  motion  of  the  electrons 
when  a  sufficiently  large  collision  frequency  is  used, 
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Plot  of  current  density,  J,  for  Case  I  in  Table  1  at  (a)  t 
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—  Schematic  of  plasma  configuration  (a)  at  early  time  (<  I  ns)  and 
(b)  at  intermediate  time  O  3  ns). 
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